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NOMENCLATURE 


area ,  cm^ 

cylindrical  surface  area  (curved  surface) ,  cm^ 
cross-sectional  area,  cm^ 

diameter  of  base  of  cone,  or  cylinder  diameter,  cm 
elastance,  elastic  load,  cmH20/liter  or  Kpa*l'^ 
cycle  frequency ,  Hz 

natural  or  resonant  cycle  frequency,  Hz 

pressure  over  atmospheric  or  ambient,  gauge  pressure, 

cmH20  or  kPa*l  ^ 

initial  air  space  height  for  water  column  container,  cm 
height  of  cone  or  pyramid,  base  to  apex,  cm 
inertance,  inertial  impedance,  cmH20/ (liter s/sec^) 
length,  axial  length  of  cylinder,  cm 
pressure,  cmH20  or  KPa 
pressure  difference,  cmH20 

flow  resistance,  resistive  impedance,  cmH20/ (liters /sec) 
specific  gravity,  dimensionless 
volume,  liters 

initial  volume  of  air  or  gas  in  water  column  container, 
liters 

volume  due  to  water  column  depression  or  movement  from 
initial  condition,  liters 
tidal  volume,  liters 


Z  »  system  impedance,  cmHaO/ (liters/ sec) 

Z*  s  complex  impedance,  cmH20/ (liters/ sec) 

|Z|  »  magnitude  of  Z*,  cmHaO/ (liters /sec) 

Zc  -  axial  length  of  cone  or  pyramid  cap,  apex  to  frustum 
plane ,  cm 

Greek  gyffit>gls 

7  =  ratio  of  specific  heats,  Cp/Cv,  dimensionless 

C  =  damping  coefficient,  dimensionless 

M  =  viscosity.  Poise 

p  =  density,  g/cm^ 

T  =  system  time  constant,  sec 

0)  =  radian  frequency,  radians /sec 
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INTRODUCTION 


Gas  flow  generated  by  a  diver  on  an  underwater  breathing 
apparatus  (UBA)  is  periodic  in  nature.  In  a  closed-circuit  UBA, 
such  as  the  Navy's  Mk-16,  breathing  impedance  comprises  flow 
resistance,  elastance  and  inertance.  These  impedances  result  in 
greater  effort  expended  to  move  gas  through  the  UBA.  From  the 
standpoint  of  scientific  understanding,  flow  resistance  is  rather 
well  understood,  but  elastance,  inertance  and  their  interaction 
are  not.  In  order  to  model  UBA  and  gain  further  fundamental 
understanding  of  elastance,  inertance  and  their  relationship  to 
UBA  breathing  impedance,  this  work  was  undertaken.  In  an  earlier 
NNRI  Technical  Report  (1)  we  studied  elastance  in  some  detail. 

In  this  work  further  issues  concerning  elastance  are  explored, 
and  inertance  is  studied  from  a  basic  experimental  viewpoint  to 
explore  ways  to  control  it. 

Elastance  and  inertance  have  important  relationships  to 
breathing  bag  design,  water  column  experiments  and  simulation  of 
UBA.  Nonlinear  elastance  may  help  to  reduce  the  elastic  load  in 
breathing  bag(s) .  The  implications  of  nonlinear  elastance  for 
impedance  testing  and  elastance  measurement  are  discussed  also  in 
this  work.  This  has  led  to  establishing  criteria  for  a  breathing 
machine . 

The  reduction  of  inertance  by  changing  the  path  of  water 
displaced  by  the  breathing  bag  also  is  explored  in  this  work,  and 
the  implications  for  natural  frequency  breathing  are  discussed. 
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ELASTANCE 


Nonlinear  ElflStange  AOd  Br.eathinq  Ms  Design 

Elastance  is  the  change  in  pressure  in  response  to  a  volume 
change  in  a  system.  Elastance  is  the  reciprocal  of  compliance 
and  has  the  units  kPa*l  \  cmH20/liter,  or  the  equivalent.  The 
pressure-volume  characteristic  of  a  closed  system  or  element 
within  that  system  can  be  obtained  by  measuring  the  pressure 
difference  when  the  volume  is  changed,  for  example,  by  a  tidal 
volume  introduced  to  the  system. 

The  pressure-volume  curve  represents  the  elastance 
characteristics  of  any  breathing  system.  When  the  elastic 
characteristics  are  linear,  the  pressure-volume  curve  is  a 
straight  line,  and  elastance  is  the  slope  of  that  line.  When 
elastic  characteristics  are  nonlinear,  the  pressure-volume 
relationship  is  curved,  and  elastance  can  be  defined  in  two  ways. 
In  the  tangent  method  elastance  is  the  slope  of  the  tangent  at  a 
point  on  the  P-V  curve,  for  example  E  =  (AP/AV)v  •  In  the  secant 
method  elastance  is  E  =  AP/AV  along  the  curve,  where  one  of  the 
P,V  points  is  always  fixed  at  a  relative  origin  (0,0).  This 
method  gives  AP/V,  which  is  the  change  in  pressure  relative  to 
the  tidal  volume  (Vj) .  We  always  use  the  secant  method  here, 
because  we  are  generally  dealing  with  tidal  volumes.  Regardless 
of  the  method  used  to  define  elastance,  the  P-V  curve  should 
always  be  obtained  for  nonlinear  systems. 
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The  elastance  of  a  water  column  is  determined  by  its 
geometry  (1).  In  a  column  with  vertical  sides  elastance  is 
essentially  linear  and  inversely  proportional  to  cross-sectional 
c^rea.  This  is  true  because  the  volume  displaced  is  equal  to  the 
area  times  the  vertical  excursion  of  the  air-water  interface. 
Since  this  vertical  excursion  is  essentially  equal  to  the 
pressure  generated  (after  appropriate  units  conversions) ,  the 
definition  of  elastance,  AP/V^  is  equal  to  (pg)h/h*A,  which 
reduces  to  (pg) /A. 

A  ncnlinear  pressure-volume  curve  can  be  generated  by  using 
a  container  with  sloping  sides,  where  cross-sectional  area  varies 
with  depth  of  the  air-water  interface.  Both  conical  and 
pyramidal  shapes  were  studied  in  the  earlier  NMRI  Technical 
Report  (1) .  These  shapes  offer  the  possibility  of  lowering 
elastic  pressures  experienced  during  breathing  by  increasing  the 
cross-sectional  area  during  a  breath.  The  equations  were 
developed  assuming  isothermal  conditions.  Adiabatic  conditions 
result  in  complications  that  make  an  analytical  solution 
impossible  for  the  cone  or  pyramid,  but  the  solution  for  the 
elastance  of  a  rigid,  fixed  volume  container  is  presented  in 
Appendix  A.  The  polytropic  index  may  prove  useful  here,  if 
desired,  as  in  Tomlinson,  et  al  (2) . 

Figure  1  illustrates  the  effects  of  orientation  and  initial 
gas  volume  on  the  elastance  of  three  conical  containers  of  equal 
size  and  geometry.  The  nonlinear  shape  of  the  P-V  curve  is 
evident  in  this  figure.  Two  cones  have  an  initial  free-space 
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Fig.  1.  Pressure-Volume  Curves  for  Conical  Bags, 
Apex  Up  and  Apex  Down 
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volume  (Vci)  corresponding  to  an  initial  free-space  height  (hj  of 
20  cm,  but  one  is  oriented  apex  up  while  the  other  has  its  apex 
down.  The  third  cone  is  oriented  apex  down  and  has  a  5  cm  free- 
space  height.  Details  of  geometry  and  nomenclature  are  given  in 
Appendix  B.  The  P-V  curves  are  all  asymmetric,  because  cross- 
sectional  area  increases  with  depth  when  the  orientation  is  apex 
up,  but  decreases  when  the  apex  is  down.  The  same  volume  change 
spread  out  over  a  larger  cross-sectional  area  gives  shallower 
depths,  hence  smaller  pressures  and  lower  elastance.  The  20  cm 
initial  free-space  height  in  the  apex  down  position  shows  highest 
elastance  in  the  positive  pressure  region  but  mid-level  elastance 
in  the  negative  pressure  region.  The  apex  up  position  with  h^  = 
20  cm  shows  low  elastance  in  the  positive  pressure  region  but 
high  elastance  in  the  negative  pressure  region,  somewhat  the 
opposite.  If  the  apex-down  cone  has  less  free-space  height,  the 
cone  is  lower  in  the  water.  For  an  initial  height  of  5  cm  in  the 
apex  down  position  pressures  are  somewhat  lower  overall  and 
elastance  is  lower,  because  the  cross-sectional  area  is 
relatively  larger  than  the  other  two  cases. 

Using  a  nonlinear  bag  in  single-bag  designs  has  no 
advantages,  because  the  reduction  in  elastance  gained  on  one 
phase  of  ventilation  would  be  lost  on  the  other.  For  example,  if 
the  apex  up  position  is  used,  elastance  is  reduced  on  exhalation, 
but  increased  on  inhalation.  The  reverse  is  true  for  the  apex 
down  position.  These  arrangements,  however,  may  be  useful  in 
physiological  testing,  for  example,  to  determine  whether 
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tolerance  for  high  elastance  is  on  the  inhalation  or  exhalation 
side,  or  at  the  beginning  or  end  of  the  tidal  volume. 

Employing  a  nonlinear  elastance  to  decrease  respiratory  load 
in  a  counter lung  can  be  demonstrated  by  a  dual  bag  system  in 
which  an  apex  down  bag  is  used  for  inhalation,  and  an  apex  up  bag 
is  used  for  exhalation.  Figure  3  shows  that  such  an  arrangement, 
when  sized  properly,  closely  approximates  the  elastance 
characteristics  of  the  Mk-15  in  the  vertical  position.  During 
inhalation  as  the  air-water  interface  of  the  inhalation  bag 
rises,  its  cross-sectional  area  increases,  thereby  reducing 
elastance.  On  exhalation,  as  the  air-water  interface  of  the 
exhalation  bag  descends,  cross-sectional  area  increases  also, 
again  reducing  elastance.  The  practical  application  of  such  a 
dual  bag  arrangement  is  complicated  by  the  necessity  of 
interconnecting  the  bags  to  provide  an  unobstructed  breathing 
circuit. 

In  summary,  to  minimize  elastance  in  a  breathing  bag,  the 
cross-sectional  area  should  be  as  large  as  practicable.  Cubic  or 
spherical  shapes  may  be  considered.  Dual  nonlinear  bag  designs 
may  be  advantageous,  but  this  will  depend  on  the  resistance  in 
the  communication  line  and  the  orientation  and  size  of  the  bags. 
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Fig.  2.  Pressure-Volume  Curves  for  Mk-15  and 
Dual  Conical  Bags 


Impedance  Testing 

Impedance  testing  with  sinusoidal  input  has  been  used  in 
physiological  studies  of  respiratory  mechanics  (3-6) .  Impedance 
is  the  ratio  of  measured  pressure  to  measured  flow  and  can  be 
represented  as  a  complex  number  with  real  and  imaginary  parts. 

Z*  =  R  +  (wl  -  E/w)j  (1) 

where  Z*  is  total  (complex)  impedance,  R  is  flow  resistance,  I  is 
inertance,  E  is  elastance  and  u  is  the  radian  frequency  (see 
Nomenclature  for  units) .  The  real  part  of  Eqn.  (1)  is  the  flow 
resistance;  the  imaginary  part  (the  reactance)  is  the  expression 
in  parentheses  multiplying  "j"  and  comprises  inertance  and 
elastance.  When  wl  =  E/«  the  reactance  is  zero,  and  the 
impedance  is  equal  to  the  flow  resistance  only.  This  defines  the 
natural  or  resonant  frequency,  which  is  given  by  =  J (E/I)  in 
radians/ sec.  The  magnitude  of  impedance  is 

|Zl  =  i/  R^  +  («I  -  E/w)^  (2) 

which  when  multiplied  by  t»ie  flow  rate  yields  the  pressure.  When 
peak  or  average  flows  are  used,  the  impedance  allows  peak  or 
average  pressures  to  be  calculated  (3) .  Equations  for  these 
components  of  impedance  in  a  variety  of  different  cases, 
including  a  water  column,  are  given  in  Appendix  C. 

The  typical  impedance  vs.  frequency  cu^ve  fox  systems  with 
linear  components  (R,  E,  I)  in  ?£•'  's  is  given  in  Fig.  3. 
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Frequency 

Fig.  3.  Impedance  Characteristics  for  System  Comprising 
Linear  Elastance,  Inertance  and  Resistance  in 
Series 
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Impedance  at  low  frequencies  is  dominated  by  elastance;  at  high 
frequencies  inertance  is  dominant.  At  the  resonant  frequency 
inertial  and  elastic  terms  cancel,  and  impedance  equals  flow 
resistance. 

The  components  of  impedance  describing  the  human  respiratory 
system  have  been  obtained  from  experimental  data  using  an 
acoustical  method  developed  by  Jackson  and  Vinegar  (6) .  In  this 
case  a  curve  similar  to  that  in  Fig.  3  is  obtained  from  pressure 
and  flow  measurements;  their  ratio  gives  the  impedance.  The 
components  of  impedance  are  extracted  by  a  multi-parameter  curve- 
fit  routine  using  a  computer.  The  technique  makes  use  of  a 
loudspeaker  and  a  controlled  frequency  sweep  in  which  the  tidal 
volume  varies  with  frequency,  having  a  larger  value  at  low 
frequencies  and  a  smaller  value  at  high  frequencies  as  a 
consequence  of  constant  power  input.  The  technique  assumes  that 
E,  I  and  R  are  independent  of  frequency. 

While  the  use  of  the  acoustic  technique  has  been  proposed  as 
a  method  of  quantifying  the  mechanical  characteristics  of  UBA 
(7) ,  the  technique  has  limitations.  When  a  UBA  is  immersed  and 
filled  to  its  operating  volume,  the  airway  pressures  can  exceed 
the  strength  of  the  speaker  cone.  Balancing  the  pressures  on  the 
other  side  of  the  cone  is  difficult  and  troublesome.  If  a 
larger,  rigid  cone  is  substituted,  excessive  force  is  required  to 
move  the  cone. 

If  elastance  is  nonlinear,  the  acoustic  technique  cannot  be 
used.  Nonlinear  elastance  is  a  function  of  tidal  volume,  which 


changes  with  frequency  during  the  test.  Knowing  the  relationship 
between  Vj  and  u  will  not  help,  because  pressure  is  the  result  of 
other  impedances  in  addition  to  elastance  and  thus  cannot  be  used 
with  Vj  to  calculate  elastance.  Obtaining  a  single  value  from 
the  impedance  curve-fit  procedure  has  no  meaning.  The  acoustical 
technique  also  uses  tidal  volumes  which  are  relatively  small 
compared  to  actual  breathing  volumes,  thus  the  elastance  measured 
by  such  a  technique  will  not  be  accurate  if  elastance  is 
nonlinear.  Even  if  the  tidal  volume  is  large  and  constant  for 
the  low  frequency  end  of  the  frequency  sweep,  the  elastance 
calculated  by  a  curve-fit  program  will  only  be  true  for  one  tidal 
volume  if  the  elastance  is  nonlinear.  Therefore,  the  test  will 
have  to  be  repeated  many  times  to  determine  the  elastance  or 
pressure-volume  characteristics  of  the  system. 

Another  difficulty  with  oscillatory  measurement  techniques 
in  general  is  that  nonlinear  elastance  will  not  be  evident  from 
the  test.  The  parameter  estimation  routine  will  curve-fit  a 
value  for  elastance  regardless  of  nonlinearity.  It  can  do  this, 
for  example,  by  using  three  points  on  the  curve  in  Fig.  3  and 
solving  Eqn.  (1)  simultaneously  for  the  three  parameters  (E,  I 
and  R)  which  are  presumed  constant.  The  "elastance"  measured  by 
this  test  may  therefore  be  totally  inaccurate. 

As  an  illustrative  example.  Fig.  4  shows  a  comparison  of 
three  curves  in  which  inertance  is  very  low  (natural  cycle 
frequency  >  1  Hz,  data  in  Appendix  D)  and  resistance  is  constant. 
The  lower  (solid)  curve  was  generated  with  linear  elastance  (5  cm 
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Impedance*  cm  water/l/e 


Frequent*  Hz 


Fig.  4.  Impedance  Characteristics  for  Linear  and 

Nonlinear  Elastance.  Curve  with  long  dashes 
indicates  test  run  with  same  Vj  —  u  dependence  as 
linear  system.  Curve  with  short  dashes  indicates 
test  run  with  different  dependence 
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H2O/ liter) .  The  middle  (long  dashed)  curve  represents  the 
nonlinear  elastance  characteristics  of  a  cone  (dimensions  in 
Appendix  D) ,  apex  down  and  Vj  =  l/u.  The  elastance  varied  from  5 
to  11.4  cm  H20/liter  at  pressures  from  1  -  30  cm  H2O, 
respectively.  The  upper  (short  dashed)  curve  represents  a  test 
of  the  same  nonlinear  system  except  that  =  2/w.  A  single 
oscillatory  test  (even  with  a  frequency  sweep)  will  not  reveal 
nonlinearity  of  elastance.  Repeating  the  test  with  different 
tidal  volume  or  tidal  volume  dependency  on  frequency  will  shift 
the  curve  as  shown.  A  curve  shift  or  change  of  impedance  at  the 
same  frequency  is  evidence  of  nonlinear  elastance,  and  a  static 
pressure- volume  test  should  then  be  performed  to  determine  the 
elastance  characteristics  from  the  P-v  curve. 

Although  the  acoustical  technique  has  too  many  difficulties 
for  application  to  UBA,  oscillatory  techniques  can  still  be  used. 
However,  any  oscillatory  technique  used  for  UBA  testing  should 
have  frequency  and  tidal  volume  independently  variable  and  have 
tidal  volumes  in  the  range  of  actual  human  breathing.  This  can  be 
done  in  a  controlled,  piston-type  breathing  machine,  such  as  that 
built  at  NMRI  for  such  studies  (7) .  In  any  oscillatory  test,  the 
test  should  be  repeated  at  least  twice  at  different  tidal  volumes 
to  check  for  potential  nonlinearity  of  elastance. 
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INERTANCE 


Water  goiumn  Experiments 

When  a  water  column  is  used  to  generate  elastance  in 
sinusoidal  breathing  experiments,  the  inertance  component  can  be 
significant  (8) .  Inertance  in  cylindrical  structures  for  various 
situations  are  compared  in  Table  1.  These  values  were  calculated 
from  equations  presented  in  Appendix  C.  Inertance  of  the  Mk-15 
(at  1  ATA  in  water)  is  seen  to  be  low,  whereas  that  of  an 
experimental  water  column  apparatus  may  be  twenty  times  higher! 
Gas  inertance  in  the  tubes  at  1  ATA  is  also  low,  but  can  increase 


TABLE  1 

Inertance  in  Cylindrical  Structures,  cmH?0/ ( 1 iters /sec^^ 


1.  Hk-15  in  water  at  1  ATA 

(D  -  30  cm,  L  =  15  cm)  .022 

2.  Cylindrical  Water  Column  Apparatus 

(A,  =  56  cm^,  L  =  25.4  cm)  .46 

3.  Air  in  a  Tube  at  1  ATA 

(A,  =  10  cm^,  L  =  180  cm,  p  =  1.2  g/liter)  .022 

4.  Air  in  a  Tube  at  Depth  (appx.  130  fsw) 

(Ax  =  10  cm^,  L  =  180  cm,  p  =  6  g/liter)  .11 


Values  computed  by 

pL/A  *  (S.G. )  (L)  /  ( .981)  (Ax)  =  cmHzO/ (liters /sec^) 
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significantly  at  depth  or  with  denser  gases  (4) ;  the  increase  in 
inertance  is  directly  proportional  to  gas  density. 

Table  2  gives  values  for  resistance,  inertance  and  elastance 
for  another  experimental  water  column  apparatus  of  larger 
diameter.  In  this  case  the  resistance  of  water  inside  the  column 
is  negligible.  As  a  result  both  inertance  and  elastance  are 
relatively  much  more  important  than  resistance.  The  natural 
cycle  frequency,  fn,  for  this  apparatus  (0.74  sec'^  or  44  breaths 
per  minute)  is  within  the  breathing  range.  The  time  constant,  t, 
is  the  inverse  of  the  radian  natural  frequency.  Breathing  on  a 
water  column  faster  than  the  natural  frequency  can  generate  peak 
pressures  in  excess  of  60  cm  H2O.  Breathing  at  the  natural 


TABLE  2 

Impedances  and  Time  Constant  for  Water  Column  Experiments 

with  Sinusoidal  Breathing 

Water  Column  Dimensions:  D  =  10  cm,  L  =  45  cm.  A,  =  78.5  cm^ 
Resistance:  R  =  S^L/A*  =  .00055  cmHjO/ (liters /sec) 

Inertance:  I  =  pL/h  =  .585  cmHaO/ (liter s/sec^) 

Elastance:  E  =  1000/A  =  12.73  cmH20/ liter 
Time  Constant:  t  =  \ (I/E)  =  .2144  sec 

Natural  Frequency:  fn  =  (27rT)  ^  =  0.74  sec'^  «  44  breaths/min 


A  conversion  factor,  981/1000,  has  been  applied  to  the  resistance 
and  inertance  formulas  to  convert  units  to  cm  H2O  and  liters. 
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frequency  minimizes  peak-to-peak  pressures  and  improves 
respiratory  comfort. 

High  inertance  may  create  difficulties  when  interpreting 
dynamic  water  column  experiments.  In  tests  with  sinusoidal 
breathing,  increases  in  the  inertial  load  cause  high  pressures 
even  though  the  elastic  load  remains  constant.  As  shown  in  Table 
1,  while  the  elastance  of  the  water  column  and  Mk-15  are  equal, 
inertance  in  the  water  column  is  much  higher  than  that  in  the  Mk- 
15  under  the  same  conditions.  Thus,  without  inertance  reduction 
the  water  column  does  not  exactly  simulate  a  diving  rig. 

Inertance  in  water  column  experiments  can  be  reduced  in 
several  ways.  One  method  is  to  shorten  the  length  of  the  tube. 

If  the  tube  is  too  short,  however,  tidal  volume  gas  will  escape. 
Placing  a  bag  at  the  end  of  a  shorter  tube  overcomes  this  problem 
and  reduces  inertance  considerably  (8).  Displacing  water 
radially  rather  than  axially  is  another  method.  This  decreases 
inertance  by  reducing  the  mass  of  water  moved  by  a  gas 
displacement  of  Vj  (the  amount  of  water  moved  is  closer  to  pV^ 
rather  than  pLA,) ,  and  second,  by  increasing  the  area  through 
which  the  mass  of  water  is  displaced  (A,  rather  than  A,)  . 

Table  3  shows  quantitatively  the  reduction  in  inertance  that 
was  achieved  by  changing  water  displacement  from  axial  to  radial. 
This  was  accomplished  in  the  laboratory  test  by  placing  holes  in 
the  sides  of  a  container  made  from  perforated  steel  sheet  bent 
into  a  conduit  with  square  cross-section  12  cm  on  a  side.  A  very 
compliant  bag  was  used  to  contain  the  gas.  A  shorter  (hence 
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lower  resistance)  flow  path  for  the  water  exists  radially  than 
exists  axially  (the  container  is  much  longer  than  it  is  wide) . 

TABLE  3 

Water  Column  Inertance  Reduction  by  Radial  Displacement 

Column  Dimensions:  D=  15  cm,  L=  60  cm,  Ajt=  177  cm^.  A,  =  2827  cm^ 

Axial  Inertance  =  pL/A,  =  mass/A,^  =  1  (60) /177  ( .  981)  = 

=  .34  cmHjO/ (liters/ sec^) 

Radial  Inertance  =  mass/Ag^  =  pV-t/ (nDL)^  =  1  (2)  (10^) /8  ( .  981)  (10®)  = 

=  .0025  cmHjO/ (liters/ sec^) 

TMs  permits  water  to  be  displaced  radially  with  a  minimum  of 
resistance.  Using  a  syringe  to  generate  pressure,  the  air-water 
interface  moved  down,  and  a  strong  current  of  water  could  be  felt 
by  placing  a  hand  near  the  sides.  Almost  nothing  was  felt  at  the 
bottom.  When  the  perforations  were  covered  with  tape  and  the 
experiment  repeated,  axial  currents  at  the  bottom  were  strong, 
and  radial  currents  were  non-existent. 

Figure  5  shows  a  pressure-time  trace  of  a  similar  test  in 
which  the  syringe  was  given  a  quick  compression  stroke  and  then 
held.  In  two  tests  the  holes  were  open,  and  in  two  others  the 
holes  were  covered  with  tape.  The  differences  were  significant. 
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TRIAL  1 


axial  (holes  taped  shut) 


radial  (holes  open) 


vertical  axis  =  pressure 
horizontal  axis  =  time 

TRIAL  2  (repeat) 


Oould  Inc..  Instrument  Systems  division 


Fig.  5.  Cessation  of  Oscillatory  Behavior  indicates 

Inertance  Reduction  in  Water  Column  Experiments 
with  Perforated-Wall  Container  (Two  Trials) 
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Lack  of  overshoot  and  oscillation  show  that  inertance  has  been 
virtually  eliminated  when  the  water  was  allowed  to  escape  through 
the  sides  rather  than  out  the  bottom. 

Brgatbinq  fiaa  Bssiqp 

Inertance  reduction  in  the  water  column  experiment  also  has 
implications  for  breathing  bag  design.  A  design  producing 
unrestricted  radial  movement  (an  expanding  sphere,  for  example) 
would  result  in  minimal  inertance.  This  occurs  because  linear 
accelerations  are  smallest  with  this  geometry.  Breathing  bags 
without  a  surrounding  housing,  e.g.  the  EX-19  design,  should 
result  in  lower  inertance  than  those  contained  within  a  housing, 
such  as  the  Mk-15.  In  the  Mk-15  holes  placed  around  the  side  of 
the  bag  housing  rather  than  at  the  back,  where  they  now  exist, 
would  promote  radial  displacement  of  water  and  reduce  inertance 
without  exposing  the  bag  to  damage.  Vfhen  water  movement  caused  by 
the  breathing  bag  is  confined  by  a  rigid  housing,  the  housing 
should  be  short  to  minimize  axial  movement. 

Tables  4  and  5  summarize  the  effects  of  axial  and  radial 
motion  using  the  Mk-15  geometry  at  pressures  of  1  ATA  and  130 
fsw.  In  the  radial  inertance  calculation,  a  2-liter  tidal  volume 
was  used.  Tidal  volume  appears  in  this  equation,  because  the 
volume  of  water  moved  is  equivalent  to  the  tidal  volume.  In  the 
axial  case,  the  volume  of  water  moved  is  governed  by  the  length 
of  the  tube  and  not  the  tidal  volume  (this  only  fixes  the 
distance  moved,  not  the  mass) .  The  elastance  of  the  Mk-15 
depends  on  orientation,  but  a  representative  value  of  4  cm 
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1 


H20/liter  was  used  (1) .  A  breathing  bag  design  promoting  the 
radial  displacement  of  water  gives  an  Inertance  more  than  an 
order  of  magnitude  lower  than  a  design  with  axial  displacement  of 
water.  At  depth  the  inertance  of  the  gas  becomes  relatively  more 
important  and  depends  on  gas  density  (depth) ,  breathing  hose 
length  and  diameter.  In  Table  5,  inertance  in  the  Mk-15  from  the 
effects  of  both  the  gas  (estimated  at  1  ATA  and  at  depth)  and  the 
water  displaced  are  shown.  Total  inertance  is  the  sum  of  both. 


TABLE  4 

Inertance  of  the  Mk-15  Breathing  Bag 
with  Axial  or  Radial  Path  at  1  ATA  in  Water 


Dimensions:  D  =  30  cm,  L  =  15  cm,  A*  =707  cm^.  A*  =  1413  cm^ 

Axial  Inertance:  pL/A,  =  1 (15) /707 ( .981)  = 

=  .0212  cmH20/ (liters /sec^) 

Radial  Inertance:  mass/A,^  =  pVj/A,^  = 

=  .0010  cmHaO/ (liter /sec^) 

Time  Constant:  t  =  J (I/E)  =  .0728  sec  (axial) 

=  .0158  sec  (radial) 

Natural  Frequency: 

fn  =  (27rT)  ^  =  2.19  sec  ^  =  131  breaths/min  (axial) 
=  10.08  sec  ^  =  605  breaths/min  (radial) 
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The  effect  of  changing  to  radial  inertance  at  depth  has  only  a 
small  ^fect  on  the  natural  frequency,  because  inertance  due  to 
gas  density  dominates  that  due  to  the  movement  of  water  arising 
from  motion  of  the  breathing  bag. 


TABLE  5 

Inertance,  Time  Constants  and  Natural  Frequencies 
in  the  Mk-15  (Simulated)  at  1  ATA  and  at  130  fsw 

Inertance .  cmHoO/  ( 1  iters /sec^l 

Water  Gas 

Axial  (baal  Radial  (baa)  at  1  ATA  at  depth 

.021  .001  .022  .110 

Time  Constant .  sec 


axial 

radial 

1  ATA 

.103 

.075 

130  fsw 

.181 

.167 

Natural  Frecruencv, 

breaths /min 

^xial 

radial 

1  ATA 

93 

127 

130  fsw 

53 

57 
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Natural  Frequency  Brga.thing 

Ideally,  the  best  UBA  design  from  a  respiratory  viewpoint:  is 
the  one  that  minimizes  peak-to-peak  mouth  pressures  for  the 
diver.  While  the  mouth  pressure  is  a  function  of  all  three 
components,  the  role  of  inertance  is  considered  to  be  minor  for 
most  practical  situations.  It  is  not  necessary,  however,  to 
eliminate  all  of  the  impedance  components  to  obtain  the  least 
encumbered  breathing.  If  the  UBA  could  be  designed  so  that  its 
natural  frequency  coincides  with  the  diver’s  breathing  frequency, 
flow  resistance  alone  contributes  to  impedance.  By  modifying 
inertance  or  elastance,  the  natural  frequency  of  a  UBA  can  be 
adjusted  to  the  breathing  frequency. 

Elastance  is  difficult  to  adjust,  but  the  possibility  exists 
that  inertance  "tuning"  can  be  useful  as  a  tool  for  reducing  peak 
pressures  and  making  breathing  more  comfortable  for  the  diver. 

As  exercise  demands  change,  the  inertance  can  be  altered  to 
"tune"  the  rig  and  lower  the  impediment  to  breathing.  For 
example,  in  the  Mk-15  inertance  can  be  reduced  by  opening  radial 
holes  instead  of  axial  holes  in  the  breathing  bag  housing,  which 
reduces  the  respiratory  load  at  high  breathing  frequency.  Table 
5  shows  some  of  the  limitations  of  this.  Controlling  inertance 
can  also  be  used  to  make  a  UBA  breathe  the  same  at  l  ATA  and  at 
depth,  compensating  for  the  increase  inertance  of  the  gas  at 
depth. 

The  effect  of  elastance  on  natural  frequency  of  a  UBA  is  the 
opposite  of  that  of  inertance.  A  high  elastance  produces  a  high 
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natural  frequency,  but  a  high  inertance  produces  a  low  natural 
frequency  (zero  inertance  produces  infinite  natural  frequency) . 
High  elastance  produces  increased  pressures  at  frequencies  below 
the  natural  frequency,  whereas  high  inertance  produces  increased 
pressures  at  frequencies  above  the  natural  frequency.  If  the 
breathing  rate  exceeds  the  natural  frequency,  higher  pressures 
are  generated  and  breathing  becomes  uncomfortable. 

The  interaction  of  elastance,  inertance  and  resistance 
should  be  a  fruitful  area  of  study  for  improved  UBA  performance 
and  diver  comfort  and  safety.  It  may  also  provide  a  basis  for 
developing  assisted-breathing  UBA's. 


SUMMARY 


Nonlinear  elastance  may  have  advantages  in  reducing  elastic 
loads  in  a  dual  breathing  bag  system  (apex  down  for  inhalation, 
apex  up  for  exhalation)  by  providing  increasing  cross  sectional 
area  on  both  phases  of  breathing. 

In  impedance  testing  nonlinear  elastance  poses  many 
problems,  and  we  recommend  static  pressure-volume  data  be 
obtained  to  determine  elastance  characteristics  of  systems  with 
unknown  or  nonlinear  elastance.  If  impedance  tests  are  used,  the 
tidal  volume  should  be  in  the  breathing  range  and  constant 
throughout  the  frequency  sweep. 

In  the  water  column  system  inertance  can  be  reduced  by 
shortening  tube  length  and  putting  a  bag  around  the  immersed  end 
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of  the  tube.  Inertance  can  be  virtually  eliminated  by  a  design 
in  which  water  is  displaced  radially  rather  than  axially,  as  the 
air-water  interface  is  moved  by  gas  displacement  inside  the 
container. 

Varying  inertance  may  be  a  useful  tool  for  adjusting  the 
natural  frequency  of  a  UBA.  By  matching  the  UBA  natural 
frequency  to  the  diver's  breathing  frequency,  the  respiratory 
load  imposed  on  a  diver  can  be  minimized,  thereby  improving  work 
capacity  and  respiratory  comfort. 
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APPENDIX  A:  DERIVATION  OF  ELASTANCE  FOR  A  RIGID  CONTAINER  UNDER 
ADIABATIC  CONDITIONS 

For  adiabatic  conditions  the  ideal  gas  law  can  be  written. 

Pi ‘Vi’  *  Pf'V/  «  constant 

The  general  equation  for  rigid  container  plus  tubing  and  syringe 
volumes  is, 

Pf  _  (V.i  +  V,  +  Veon)’ 

(V.i  +  V,  Veo„  -  Vt)" 

or 

1  +  h/Pi  =  -  Vt)’ 

and 

V.«i"  -  (V,«i  -  Vi)’’ 

h  =  Pi  •  _ 

(V.„i  -  Vt)" 

Taking  the  derivative  with  respect  to  Vj, 

d  V.„i’ 

dh/dVi  =  E  =  Pi  • _  _  -  1 

dV,  (V.„i  -  Vt)’ 

=  Pi-V.„J(-7(V.«i  -  Vt)  ■■'*'( -1) 

=  7PiV.u»iV  (V,u„.i  - 

for  V,„i>  Vx, 

dh/dVx  =  7Pi/V.„„.i  =  7  Po/Vo- 
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APPENDIX  B;  CONE  DIMENSIONS  AND  DATA  FOR  FIGURES  1  AND  2. 

The  equations  for  a  conical  shape  are: 

Cone ,  Apex  Up 

Vj  =  (itdVish^)  (h/1000)  [h^  +  3h{hi  +  2c)  +  3(hi  +  Zc)^]  (1) 

Cone ,  Apex  Down 

Vt  =  (ttdViSH^)  (h/1000)  [h^  -  3h(H  -  hi)  +  3  (H  -  hi)^]  (2) 

Where  V,  is  tidal  volume,  H  is  height  dimension  of  the  cone  from 
base  to  apex,  h  is  the  depth  {+)  or  height  (-)  of  water  in  the 
container  relative  to  that  outside  the  container  (h  is  also  gauge 
pressure),  hi  is  the  initial  height  of  gas  spcce  in  the 
container,  is  the  height  of  the  top  (apex)  part  of  the  cone 
removed  to  create  a  frustum,  and  D  is  the  cone  base  diameter. 
These  dimensions  are  illustrated  in  Fig.  6  for  both  cases. 
Equations  1  and  2  are  dimensional  with  Vj  in  liters  and  all 
distances  in  cm. 

Cone  Dimensions:  D  =  20  cm,  H  =  50  cm,  z^  =  15  cm 


28 


EX  DOWN 


APPENDIX  B  (cont.) 
A.  Figure  i 
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3.  hi  20  cm,  cone  apex  up 
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B.  Figure  2 

Data  for  Mk-15  from  reference  (1) 

Data  for  composite  curve:  negative  pressures  from  Fig.  1, 
hi  =  5  cm,  cone  apex  down;  positive  pressures  from  Fig.  1, 
hi  =  20  cm,  cone  apex  up 
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APPENDIX  C:  SUMMARY  OF  EQUATIONS  FOR  IMPEDANCE  COMPONENTS 

Peslln  and  Fredberg  (3)  sxuninarize  basic  impedance 
relationships  for  sinusoidal  inputs.  The  impedance  (Z)  is  a 
complex  quantity  with  the  real  part  being  the  resistance  and  the 
imaginary  part  combining  inertance  and  elastance  according  to 
Eqn.  3.  In  the  following,  t  is  the  system  time  constant,  fn  is 
the  natural  or  resonant  frequency,  and  ^  is  the  damping 
coefficient.  If  the  damping  coefficient  is  1.0,  the  system  is 
said  to  be  critically  damped  (no  oscillations  occur  in  response 
to  a  system  disturbance) .  If  the  damping  coefficient  is  less 
than  1,  the  system  is  underdamped,  and  oscillatory  behavior  is 
expected.  If  the  damping  coefficient  is  greater  than  1,  the 
system  is  overdamped,  and  a  smooth  time  response  rather  than 
oscillatory  behavior  is  expected. 

Resistance:  R  =  f (geometry,  flow,  w,  p,  p) 

Elastance:  E  =  7P0/V0  (adiabatic,  rigid  container) 

=  Pi/Vi  (isothermal,  rigid  container) 

Inertance:  I  =  npL/A  (n=l  for  flat  velocity  profile, 

n=4/3  for  parabolic  profile) 

I  =  mass/A^ 

Time  Constant:  t  =  \ (I/E) 

Natural  Freouency:  f„  =  (2itt)'^ 

Damping  Coefficient:  2fT  =  R/E,  Thus  f  =  R/2j (IE) 
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For  a  Cylindrical  Water  Column  with  Nonporous  Sides: 

R  »  S/^L/A*  (for  laminar  flow,  linear  resistance) 

E  =  1/A 

I  =  pL/A  =  mass/A^ 

Units  conversion  factors  must  additionally  be  used  to  obtain  the 
units  of  cm  HaO  and  liters. 

In  water  column  experiments,  resistance  and  inertance  cannot 
be  defined  exactly.  Both  are  a  function  of  mass  which  is  not 
constant,  since  the  level  of  water  inside  the  cylinder  changes 
with  pressure. 


I 

i 


33 


APPENDIX  D:  IMPEDANCE  TEST  SIMULATION  AND  DATA  FOR  FIGURE  4 


1.  Data  for  Figure  4:  Linear  and  Nonlinear  Elastance  in 

Impedance  Tests 
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Cone  with 

hi  =  10 

cm,  D 

=  20  cm. 

H  =  50 

cm. 

2c  =  15 

cm 

Vx  =  4.19  X  10'^  (h^  -  120  +  4800  h) 


R  =  3 

Z  =  /  (9  +  eV«^) 

Zi  is  linear,  Zz  and  Zj  are  nonlinear  (same  model,  different 
frequency-tidal  volume  dependence) . 
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